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We present a measurement of the top quark mass and of the top-antitop (tt) pair production 
cross section using pp data collected with the CDF II detector at the Tevatron Collider at the Fermi 
National Accelerator Laboratory and corresponding to an integrated luminosity of 2.9 fb _1 . We 
select events with six or more jets satisfying a number of kinematical requirements imposed by 
means of a neural-network algorithm. At least one of these jets must originate from a b quark, 
as identified by the reconstruction of a secondary vertex inside the jet. The mass measurement 
is based on a likelihood fit incorporating reconstructed mass distributions representative of signal 
and background, where the absolute jet energy scale (JES) is measured simultaneously with the 
top quark mass. The measurement yields a value of 174.8 ± 2.4(stat+JES) j^j^syst) GeV/c 2 , 
where the uncertainty from the absolute jet energy scale is evaluated together with the statistical 
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uncertainty. The procedure also measures the amount of signal from which we derive a cross section, 
o t t = 7.2 ± 0.5(stat) ± l.O(syst) ± 0.4(lum) pb, for the measured values of top quark mass and JES. 



PACS numbers: 14.65.Ha, 13.85.Ni,13.85.Qk 

I. INTRODUCTION 

Since its early measurements, the large value of the 
top quark mass (Mt op ) has represented a really strik- 
ing property of this particle, giving to the top quark a 
special position within the standard model (SM) and 
suggesting also possible links to new physics |l| . In 
fact, apart from being itself a fundamental parameter 
of the SM, M top is by far the largest mass among the 
ones of the observed fermions, and this makes the top 
quark contribution dominant in higher order correc- 
tions to many observables. Therefore M top plays a 
central role in checking the consistency of theoretical 
predictions of the SM. The higher order corrections 
apply also to the W boson propagator, and therefore 
affect the calculated value of the W mass, Mw- As 
the latter depends logarithmically on the mass of the 
Higgs boson, precise measurements of M\y and M top 
allow setting indirect constraints on the value of the 
mass of this fundamental, but still unobserved parti- 
cle 0. Moreover, possible contributions due to some 
unknown physics might also be constrained. Finally, 
the present value of M top makes the Yukawa coupling 
to the Higgs field of 0(1) and this could indicate a 
special role of the top quark in the mechanism of elec- 
troweak symmetry breaking. 

All these reasons make the accurate knowledge of 
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M top a really important issue, but the same is true 
for the measurement of the ti production cross section 
(<J t t), both as a test for physics contributions beyond 
the SM and as a test of current next-to-leading-ordcr 
(NLO) QCD calculations 0. Usually, measurements 
of a t t rely upon event counting and are performed as- 
suming an a priori value for M top - The technique 
used here allows the simultaneous measurement of 
both these important and related properties of the 
top quark. 

At the Tevatron Collider at Fermi National Accel- 
erator Laboratory, top quarks are produced mostly in 
pairs. In the SM the top quark decays into a W boson 
and a b quark almost 100% of the time, and the topol- 
ogy of the final state resulting from a ti event depends 
on the hadronic or leptonic decay of the two final-state 
W bosons. In this paper, we consider events char- 
acterized by a multijet topology (all-hadronic mode) 
with no energetic leptons. This ti final state has the 
advantage of a large branching ratio (« 4/9) and of 
having no undetectable final-state particles. The ma- 
jor challenge of this channel is the large background 
from QCD multijet production, which dominates the 
signal by 3 orders of magnitude after the application 
of a specific online event selection (trigger). To in- 
crease the purity of the candidate sample, require- 
ments based on the kinematical and topological char- 
acteristics of SM ti events are expressed in terms of 
an artificial neural network and applied to the data. 
Further improvement is then obtained from the re- 
quirement of at least one jet identified as originating 
from a b quark using a secondary vertex &-tagging al- 
gorithm. Simulations predict that a clear ti signal will 
thus become visible over background in the selected 
data sample, and the measurement of the top quark 
mass and the ti cross section is made possible in spite 
of the overwhelming QCD multijet production. 

A reconstructed top quark mass is determined by 
fitting the kinematics of the six leading jets in the 
event to a ti final state. This variable, denoted as 
mj ec , does not strictly represent a measurement of 
M top , but its distribution obtained by a sample of ti 
events is sensitive to Mt op itself. The jet energy scale 
(JES) is a factor representing the set of corrections 
needed to obtain a better estimate of the energy of a 
parton starting from a jet reconstructed by clusters in 
the calorimeter. The default JES used in simulated 
events is obtained by a tuning to the data, but pos- 
sible discrepancies between data and simulation lead 
to an uncertainty on this value. The strong correla- 
tion existing between the m\ ec distribution and the 
JES implies therefore a corresponding uncertainty on 
M t0 p- However, the JES can be calibrated using the 
selected samples of ti candidate events, where a second 
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variable, wi^> c , is reconstructed by the four- momenta 
of the jets assigned to the W bosons. This variable 
is related to the well-known value of the VF-boson 
mass, and the JES can be adjusted in such a way that 



both the 



and the : 



distributions for simulated 



events match the observed data. The inclusion of this 
procedure, usually referred to as in situ calibration, 
enables a significant reduction of the systematic un- 
certainty associated with the inaccurate knowledge of 
the JES, and represents an important improvement of 
the work described in this paper with respect to the 
previous CDF analysis by a similar method [1] . 

The m\ ec and rriyy distributions are reconstructed 
in two separate samples of selected data events, de- 
fined by the presence of exactly one and two or more 
6-tagged jets respectively. The data are then com- 
pared to corresponding distributions expected from 
background and tt events simulated with various val- 
ues of the top quark mass and of the JES to fit for 
these parameters. In addition, the fitted signal yields 
are used to derive a measurement of the tt production 
cross section. 

The results reported here are based on data taken 
between March 2002 and April 2008, corresponding 
to an integrated luminosity of 2.9 fb^ 1 . This mea- 
surement complements other recent determinations of 
the top quark mass and tt cross section by CDF and 
DO [5], Q in other final states, and improves the latest 
CDF measurements in the same channel (3, 0] • 

The organization of the paper is as follows: Sec- 
tionlnlcontains a brief description of the CDF II detec- 
tor. The trigger and the neural-network-based sample 
selection are discussed in Sec. lIIIl along with the iden- 
tification of jets initiated by b quarks (b jets). Sec- 
tions lIVI and|V] present the simulated signal samples 
and the data-driven method we use for estimating the 
background from multijet data. Section lVTl describes 
how the fundamental variables m\ cc and m^y are re- 
constructed, while in Sec. lVIII we present the final re- 
quirements to define the samples of events used in the 
measurement. The parametrization of the dependence 
of the distributions of reconstructed variables on the 
values of the top quark mass and the jet energy scale 
are described in Sec. lVIITAl The fit to the experi- 
mental distributions and its calibration are described 
in Secs. lVIII Bl and llXl respectively. SectionlXl details 
the study of the systematic uncertainties on the mass 
measurement, that is then reported in Sec. lXIl We 
describe in Sec. lXIII the measurement of the tt cross 
section. 



II. THE CDF II DETECTOR 

The CDF II detector [8] is an azimuthally and 
forward-backward symmetric apparatus designed to 
study pp collisions at the Tevatron. A cylindrical 
coordinate system is used where 9 is the polar an- 
gle to the proton beam direction at the event vertex, 



(f> is the azimuthal angle about the beam axis, and 
pseudorapidity is defined as 77 = — In [tan(0/2)]. We 
define transverse energy as Et — E sin 9 and trans- 
verse momentum as pr — psin(9 where E is the en- 
ergy measured by calorimeters, and p is the magni- 
tude of the momentum measured by a tracking sys- 
tem. The detector consists of a magnetic spectrome- 
ter surrounded by calorimeters and muon chambers. 
The charged particle tracking system is immersed in 
a 1.4 T solenoidal magnetic field with axis parallel 
to the beamline. A set of silicon microstrip detectors 
provides charged particle tracking in the radial range 
from 1.5 to 28 cm, while a 3.1 m long open-cell drift 
chamber, the central outer tracker (COT), covers the 
radial range from 40 to 137 cm. In combination the sil- 
icon and COT detectors provide excellent tracking up 
to about pseudorapidities \n\ < 1.1, and with decreas- 
ing precision up to < 2.0. Segmented electromag- 
netic and hadronic calorimeters surround the tracking 
system, and measure the energy deposit of particles 
interacting in the calorimeters. The electromagnetic 
and hadronic calorimeters are lead-scintillator and 
iron-scintillator sampling devices, respectively, cover- 
ing the range \r)\ < 3.6. They are segmented in the 
central region (|?7| < 1.1) in towers of 15° in azimuth 
and 0.1 in ij, and the forward region (1.1 < \rj\ < 3.6) 
in towers of 7.5° for \rj\ < 2.11 and 15° for |r?| > 2.11, 
while the coverage in 1 77 1 increases graduall y fr om 0.1 
to 0.6 . The electromagnetic calorimeters [9lTTo| are in- 
strumented with proportional chambers (at large an- 
gles) or scintillating strip detectors (in the forward 
regions) , which measure the transverse profile of elec- 
tromagnetic showers at a depth corresponding to the 
expected shower maxima. Drift chambers located out- 
side the central hadronic calorimeters and behind a 
60 cm iron shield detect muons with \rj\ < 0.6 11 1. 



Additional drift chambers and scintillation counters 
detect muons in the region 0.6 < \rj\ < 1.5. Mul- 
ticell gas Cherenkov counters [l2| with a coverage of 
3.7 < < 4.7 measure the average number of inelas- 
tic pp collisions and thereby are used to determine the 
luminosity. 



III. MULTIJET EVENT SELECTION AND b 
TAGGING 



The final state of all-hadronic tt events is charac- 
terized by the presence of at least six jets from the 
decay of the two top quarks, where additional jets 
might come from initial- or final-state radiation (ISR 
or FSR). Events having such a topology are collected 
using a multijet trigger which relies on calorimeter 
information. Subsequently, jets are identified during 
event reconstruction by grouping clusters of energy in 
the calorimeter using a fixed-cone algorithm with a 
radius of 0.4 in r\-(\> space [l3j]. After a preliminary se- 
lection of multijet events, a neural-network selection 
based on relevant kincmatical variables is used to fur- 



G 



ther improve the purity of the sample. 



A. Multijet Trigger 

The CDF trigger system has three levels. The first 
two levels consist of special-purpose electronic circuits 
and the third one of conventional programmable digi- 
tal processors. At level 1, the trigger requires the pres- 
ence of at least one calorimeter tower with transverse 
energy Eft w > 10 GeV. At level 2 the total trans- 
verse energy, obtained as the sum over all calorimeter 
towers, J2 E T W y must be ^ 175 GeV. Moreover, the 
presence of least four clusters of towers, each with 
transverse energy E T lus > 15 GeV, is required. Fi- 
nally, the third trigger level confirms the level 2 se- 
lection using a more accurate determination of the jet 
energy, requiring four or more reconstructed jets with 
Et > 10 GeV. Approximately 14 x 10 6 events satisfy 
the trigger requirements, corresponding to an events 
signal-over-background ratio (S/B) of about 1/1200, 
assuming a theoretical cross section of 6.7 pb for a top 
quark mass of 175 GeV/c 2 |3|. 



B. Preselection and Topology Requirements 

Events satisfying the trigger requirements are re- 
constructed in terms of their final-state observables 
(tracks, vertices, charged leptons, and jets). We retain 
only those events that are well contained in detector 
acceptance, requiring the primary event vertex [7JJ to 
lie inside the luminous region (\z\ < 60 cm). We 
remove events having well identified energetic elec- 
trons or muons as defined in [l5j , namely electrons 
with E T > 20 GeV and muons with p T > 20GeV/c. 

In order to have jets matching as accurately as pos- 
sible to the hard scattering partons, we correct jet 
energies for detector response and multiple interac- 
tions 16]. First, we consider the 77-dependence of de- 
tector response and energy loss in the uninstrumented 
regions. Then, after accounting for the small extra 
energy deposited by multiple collisions in the same 
beam-beam bunch crossing, a correction for calorime- 
ter nonlinearity is applied so that the jet energies are 
equal, on average, to the energy of the particles in- 
cident on the jet cone. The total uncertainty on the 
estimate of the original parton energy, where all un- 
certainties for the individual corrections are added in 
quadrature, varies from 8% to 3% with jet transverse 
energy increasing from 15 GeV to 50 GeV, and re- 
mains approximately constant at 3% above 50 GeV. 
Jets with \r/\ < 2 and Et > 15 GeV, after all correc- 
tions are applied, are selected for further analysis. 

As the uncertainty on the missing transverse en- 
ergy, $t [13 1 increases proportionally to y/^Er , 

its significance is defined as , T , where the Et is 

corrected for any identified muons, while ^ Et is ob- 



tained by summing the Et's of all the selected jets. 

E i 
We then require that , T be < 3 GeV 2 to select 

\JY.Et 

events with small $t- At this stage, called preselec- 
tion, we are left with about 8.2 x 10 6 events. 

As the topology of the candidate events is de- 
termined by the jet multiplicity, we define the sig- 
nal region by selecting events with a number of jets 
6 < iVjcts < 8 and we also require jet pairs to be sepa- 
rated by at least 0.5 units in the rj-(f> space. The num- 
ber of events passing these additional requirements is 
1.671 x 10 6 with an expected S/B of approximately 
1/430. 



C. Neural-Network-based Kinematical Selection 

To further improve the purity of the signal sam- 
ple, we use a multivariate approach and take advan- 
tage of the distinctive features of signal and back- 
ground events through a neural network, which takes 
into account the correlations between the kinemati- 
cal variables which enter as input nodes in the net- 
work. The network uses the mlpfit package [TJI as 
implemented by ROOT [20| through the TMultiLayer- 
Perceptron class. 

A first set of 11 global variables, summarized in 
Table[H have already been proven to be effective Q in 
reducing the QCD background. Studies performed for 
this analysis on the jet development in the calorime- 
ter have indicated that a good discrimination between 
quark-initiated and gluon-initiated jets can be accom- 
plished with 77 moments (M v ) and <f> moments (Ma,) 
of a jet, which are defined as 



M n = 



\ 



^ — ' Et 

tow 



"^7tow 



(1) 



and 



Ma, 



U 



iT'tow 

,2 

P fto 
HiT 



(2) 



where Et, rj, and <p are j respectively, the transverse 
energy, the pseudorapidity and the azimuthal angle of 
the jet, while E^f™ is the transverse energy deposited 
in the calorimeter towers belonging to the jet. 

We remove possible biases coming from Et distri- 
butions, which might differ in signal and background 
events, by deconvoluting the Et dependence through 
a rescaling of all moments to a common reference value 
of Et = 50 GeV. We obtain what we call scaled mo- 
ments : 



M° = M n 



/?(50 GeV) 



(3) 
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and 



Ml = M*x 



#(50 GeV) 



(4) 



where f^(Er) and /^(£r) are the functions that fit 
the profiles of M„ vs Et and of M^ vs Ex in quark- 
initiated jets from simulated ti events. 

These scaled moments are quite different for jets 
coming from a quark or a gluon in simulated ti events. 
Such a behavior has been verified in data events where 
the jet origin is well known. To take advantage of 
the large number of jets in a ti event, we consider 
the geometric average of the rj moments and of the cj> 
moments, see Fig.[TJ evaluated using all jets which are 
not identified as coming from a heavy quark by the 
criteria explained in Sec. MI Dl 

The 13 variables are used as inputs to a neural net- 
work with two hidden layers with 20 and 10 hidden 
nodes, respectively, and one output node. The net- 
work is trained on same-size samples of signal and 
background events with 6 < iVjets < 8 (about half a 
million events). In order to model the signal we use 
the PYTHiAv6.2 leading-order (LO) Monte Carlo 
generator with parton showering followed by a simula- 
tion of the CDF II detector. The reference top quark 
mass chosen for the training is M top — 175GeV/c 2 . 
The background is obtained from the multijet data 
events themselves, since the signal fraction is expected 
to be very small before applying the neural-network 
selection. The value of the output node, N out , is 
the quantity we use as a discriminator between sig- 
nal and background, and is shown in Fig. [2] for the 
6 < IVjets < 8 sample. 



TABLE I: Input variables to the neural network. 



Variable 



Description 



C 
A 

M max 
M max 

e: ' 

(E T ) 
(M°) 
(Ml) 



T 
*.2 



Scalar sum of selected jets Et 

As above, except the two highest--&r jets 

Centrality 

Aplanarity 

Minimum dijet invariant mass 
Maximum dijet invariant mass 
Minimum trijet invariant mass 
Maximum trijet invariant mass 
Et sin 2 6* for the highest-Br jet 
Et sin 2 9* for the next-to-highest- Et jet 
Geometric mean over the remaining jets 
Geometric mean over the untagged jets 
Geometric mean over the untagged jets 




0.12 0.14 0.16 

Average^ moment 



FIG. 1: Geometric average of the r\ scaled moments 
({M%), upper plot) and of the <f> scaled moments ((Ml), 
lower plot) for QCD multijet (solid histogram) and simu- 
lated it (dashed histogram) events with 6 < A/jets < 8. 




FIG. 2: Neural-network output iV ou t for QCD multi- 
jet (solid histogram) and simulated ti (dashed histogram) 
events with 6 < Nj ets < 8. Histograms are normalized to 
unity. The neural-network implementation that we use in 
the TMultiLayer-Perceptron produces an output which is 
not strictly bound between and 1. 



D. Tagging b quarks 

In order to enrich the ti content in the event 
sample, we use a 6-tagging algorithm based on sec- 
ondary vertex reconstruction as described in detail 
in [3, [22j ■ The algorithm identifies a jet likely to 
contain a hadron with a b quark by reconstructing 
its decay vertex with at least two high-quality tracks 
with hits in the silicon vertex detector. A 6-tagged 
jet {tag, in brief) must have an associated secondary 
vertex with a displacement from the primary ver- 
tex in the transverse plane larger than 7.5 times the 
transverse-displacement resolution. This is evaluated 
for each secondary vertex, but its typical value is 
about 190 fim. The tagging efficiencies for jets coming 
from the fragmentation of b or c quarks are corrected 
in simulated events according to the efficiency seen in 
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the data, by a factor 0.95 ± 0.04, both for b jets and c 
jets. These factors are described in detail in [14J. 



IV. EVENT SIMULATION 

The standard model tt events used to study the 
event selection and to check the performance of the 
method (Sec.lTX|) are simulated using pythia v6.2 |2lj . 
Samples generated with input values of the top quark 
mass, MJo p , ranging from 160 to 190 GeV/c 2 are con- 
sidered and, for each sample, the event selection is re- 
peated by varying the JES from its default value [16[ . 
The displacement, denoted as AJES, is measured rel- 
ative to the uncertainty, ctjes, on the default value 
itself, so that the value of JES applied to simulated 
events is increased by AJES • ctjes with respect to 
the default. To test the method, input values AJES m 
ranging from —3 to +3 are considered. 

Different generators and different values for the 
model parameters are used to estimate the systematic 
uncertainties, as described in Sec. IXl 



V. BACKGROUND ESTIMATE 

The background for the tt multijet final state comes 
mainly from QCD production of heavy-flavor quark 
pairs (bb and cc) and from false tags of light-flavor 
quark jets. Other standard model processes such as 
W/Z+jets have a smaller production cross section and 
small acceptance due to the selection requirements. 

Given the large theoretical uncertainties on the 
QCD multijet production cross section, a more accu- 
rate background estimate is obtained from the data, 
rather than from Monte Carlo simulations. A tag 
rate per jet, defined as the probability of tagging a 
jet whose tracks are reconstructed in the vertex de- 
tector (fiducial jet), is then evaluated in a sample of 
events with exactly four jets passing the preselection 
and therefore still dominated by the background (S/B 
~ 1/5000). The rate is parametrized in terms of vari- 
ables sensitive to both the tagging efficiency for heavy- 
flavored objects and the probability of false tags : the 
jet Et, the number of tracks reconstructed in the sili- 
con vertex detector and associated with the jet, N^ k , 
and the number of primary vertices reconstructed in 
the event, N VCIt [4]. By definition, the tag rate es- 
timates the probability that a fiducial jet having, on 
average, the characteristics of jets from background 
events is tagged. Its average value is about 3.7%, with 
negligible uncertainty. However, direct exploitation 
of the tag rate to predict the number of background 
events with exactly a given number of tags would give 
incorrect numbers. This happens because, by con- 
struction, this rate is the ratio between the number of 
tagged jets and the number of fiducial jets in a whole 
sample of events. Possible correlations among jets in 
the same event are not considered. As heavy-flavor 



quarks come in pairs in QCD background, the proba- 
bility to tag a pair of jets in the same event is therefore 
larger than the simple product of the tag probabilities 
of individual jets given by the tag rate. 

To account for this we introduce correction factors 
to obtain a better estimate for the number of 1-tag and 

> 2-tag background events. These factors are derived 
in a control sample dominated by the background 
(events with six, seven or eight jets and A ou t < 0.25, 
with S/Bw 1/1300 for one tag and S/Bw 1/400 for 

> 2 tags) as the ratio between the observed number 
of events with n tags (with n = 1, 2, 3) and the av- 
erage expectation obtained by using the tag rate to 
evaluate the probability for each event to have the 
same number, n, of tagged jets. These factors rep- 
resent, therefore, average corrections to the probabil- 
ity for a possible tag configuration, that is for the as- 
sumption that among the fiducial jets in an event of 
the sample selected before the 6-tagging requirements 
(pretag sample) only a given subset is actually tagged 
when the algorithm is applied. Their average values 
are 0.94, 1.48, and 2.46 for events with one, two, and 
three tagged jets, with relative statistical uncertain- 
ties of 0.4%, 1.1%, and 5.1% respectively. Similarly 
to the tag rate, these corrections should be valid for 
events with the characteristics of background events. 

The accuracy of our modeling of the background 
processes is verified in control samples, i.e. on events 
with higher values of N ou t and therefore with a larger 
fraction of signal events and with possible different 
kinematics and background composition. As the back- 
ground prediction is performed using the data in the 
pretag sample, the presence of tt events must also be 
taken into account. Therefore a correction is applied 



to derive a better evaluation, n 



(b, exp) ■ 



of the back- 



ground normalization from the raw estimate nr^ raw ) 
directly obtained by the corrected tag rate matrix. 
This correction must subtract the contribution, n!?, 
coming from applying the matrix to signal events and 
included in n.({,, raw) • Denoting by A Q b s the number of 
events observed in the data sample, by n t j the number 
of signal events in this sample and assuming that the 
excess of events with respect to the expected back- 
ground is totally due to the signal, the correction can 
be written as : 

tr 

""(6, exp) raw) 

tr 

— n (b, raw) ' n tt 

n-tt 
n tr - 

= ^(6, raw) — • (^obs - «-(fc,cxp)) > (5) 

which, with R t i = n^/nu, gives 

n (b, raw) - Rtt ■ N obs , . 

"(b,cxp) = \^R~ ■ W 

R t t can be inferred from simulated events and 
amounts to 0.314 ± 0.003 (stat) [0.067 ± 0.0014 (stat)] 
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for 1-tag (> 2-tag) events. Further possible discrep- 
ancies between the observed and expected number of 
events are considered as due to the modeling of the 
background and accounted for as a systematic uncer- 
tainty. 



VI. MASS RECONSTRUCTION 

The simultaneous measurement of the top quark 
mass and the JES is based on the reconstruction, event 
by event, of both the top quark and the W masses 
through a constrained fitting technique. The shapes of 
the distributions obtained by this procedure are sen- 
sitive to the values of both M top and JES. Therefore, 
for simulated events, they are built using samples cor- 
responding to the different input values of Af top and 
AJES listed in SeaEl 

Moreover, given the different resolution in the re- 
constructed top quark mass and the W boson mass, 
and also the different S/B which can be achieved by 
requiring events with exactly one or > 2 tags, two sets 
of distributions are separately derived in these sam- 
ples. 



A. Reconstructed top quark mass 



reconstructed minimizing the following \ 2 function : 
2 ( m \V ~ M w i 



(1) 



(m {2 }-M w ) 2 



33 



r 2 

1 w 

I (!) 

\ m jjb ~ m 



r 2 

1 w 



) 2 , ( 



,(2) 



m jjb ~ m 



+ 



r 2 



\Pr,i Pr,i 



r 2 



i=i 



(ii) 



The minimization procedure is performed with respect 
to seven parameters, i.e., the reconstructed top quark 
mass m\ cc and the transverse momenta of the 
six jets, which are constrained to the measured value 
p™ as within their known resolution er^. The invari- 
ant masses of the jet doublets assigned to light-flavor 



quarks coming from a W, 



,(1,2) 



and of the trijet 



systems including one doublet and one of the jets as- 

(1 2) 

signed to b quarks, JtoLj, , are evaluated by the trial 
momenta of jets at each step of the minimization. On 
the contrary, the measured mass Mw and the natu- 
ral width T\v of the W boson as well as the assumed 
natural width of the top quark, T t , are kept constant 
to 80.4GeV/c 2 , 2.1 GeV/c 2 and 1.5GeV/c 2 respec- 
tively m ni . 

The permutation of jets which gives the lowest % 2 
value is selected, and the corresponding fitted value of 
m\ ec enters an invariant mass distribution {template) 
which will be used for the M top measurement. 



For each event we 



determine a reconstructed 
top quark mass, m\ ec 1 from the four- momenta of 
selected jets. Sixteen equations can be consid- 
ered to connect the four-momenta of the two top 
quarks and their decay products according to the 
tt—> bbW + W~ — !• W>qi<Z2<Z3<Z4 hypothesis: 



Pt 

a 

Pi 

a 

Pw+ 



a 

Pw+ 
Pw- 



'91 



p. 



w- 



P' 



■p%1 
■pt 



(7) 
(8) 
(9) 
(10) 



with fj, = 0, 1, 2, 3. There are 13 unknown quanti- 
ties, i.e., the unknown top quark mass and the three- 
momenta of the top quarks and of the W bosons, so 
the kinematics of the events are overconstrained. 

The fit is performed using only the six highest- Et 
jets (leading jets) of the event and considering their 
possible assignments to quarks of a ti final state. The 
total number of different permutations giving two dou- 
blets of jets corresponding to the W bosons and two 
triplets of jets corresponding to the top quarks is 90. 
Since we require the presence of b tags, assigning the 
tagged jets only to b quarks reduces this number to 
30 for 1-tag events and 6 in case of two or more b 
tags . 

For each permutation the kinematics of the event is 



Reconstructed W mass 



Reconstructing the mass of W bosons by using dijet 
systems represents a possibility to obtain a variable, in 
principle, insensitive to M top which allows, therefore, 
an independent determination of JES. 

To build the m^> c distributions we use the same 
procedure and \ 2 expression considered for m r t ec , but 
now the W-boson mass is also left as a free parameter 
in the fit (i.e. Mw becomes m^ c ). Again, for each 
event, the value of rn^ corresponding to the permu- 
tation of the jet-to-parton assignments with the lowest 
X 2 enters the distribution. 

Using different fits in the reconstruction of mj ec 
and 



MI 



can lead to selecting different assignments 
of jets to partons for the two variables in the same 
event. This is not a problem as the same procedure is 
followed both on data and simulated events. Recon- 
structing the top quark mass using a constant value of 
Mw, as decribed in Sec. I VI Al improves the resolution 
of the distributions and therefore the determination of 
the true value of M top . The correlations between the 



values of 



and 



Mr 



in the same event are taken 



into account in the calibration of the likelihood fit used 
for the measurement fSec lVHiBl) . 
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C. Background templates 



In order to reconstruct data-driven background 
templates we apply the kinematical fitter to the sam- 
ple of events passing the neural-network selection, but 
before the requirement of tagged b jets. 

The same procedures described in Sees. I VI Al 
and IVIBI are repeated on these events assigning fidu- 
cial jets to b quarks and then looping over all possi- 
ble assignments of other jets to the remaining quarks, 
performing the fit for each permutation and selecting 
the reconstructed m\ cc and m^ c values correspond- 
ing to the best x 2 ■ These values then enter the tem- 
plates weighted by the corrected probability of the as- 
sumed tag configuration; see Sec. |V] As for the nor- 
malization, the background distributions also need to 
be corrected for the presence of signal in the pretag 
sample by subtracting the contribution from tt events. 
The shape of this contribution is obtained from simu- 
lated samples and depends on the assumed M top and 
JES, while the normalization is given by the difference 
n (b, raw) — n (&,oxp): as described in Sec.lVl 

In order to check how well our modeling describes 
the background, we consider events in control regions 
defined by the N ou t value, in ranges where the signal 
presence after tagging is still very low. In these regions 
the templates, i.e. the main elements of our measure- 
ment, are reconstructed by the procedure described 
in the previous sections, both for the signal and the 
background, as well as other important distributions 
like iV ou t and the \ 2 01 the fit used to build the m\ ec 
templates. These distributions are then compared to 
observed data, taking into account the contribution 
from signal events. The agreement is generally good 
in all the control regions, and this confirms the relia- 
bility of the background model. 

Figures [3] and 0] show, as examples, distributions 
of m\ cc and m^ c in one of the control regions for 
1-tag and > 2-tag events, where the sum of signal 
and background is compared to the same distribu- 
tions reconstructed in the data. In these plots the 
integral of the signal distributions corresponding to 
M top = 175GeV/c 2 and the default value AJES=0 
have been normalized to the difference between the 
observed data and the corrected expected background. 
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FIG. 3: Histograms of the reconstructed top quark mass 
mj ec for 1-tag events, upper plot, and > 2-tag events, 
lower plot, are shown in a control region defined by 
0.75 < Nout < 0.85. Along with the data are plotted 
the expected background and the signal contribution for 
M to p = 175GeV/c 2 and the default value AJES = 0, nor- 
malized to the difference between the data and the back- 
ground. The value of the purely statistical \ 2 probability 
is reported on each plot. 



VII. EVENT SAMPLES 



In order to obtain the best performance from 
our method, we performed sets of pseudoexperiments 
(PEs) to find out which requirements on the values of 
7V out and of the x 2 used to obtain the m\ cc values min- 
imize the statistical uncertainty on the top quark mass 
measurement. The procedure is similar to the one out- 
lined in Sec. lIXl with a binned version of the same like- 
lihood. It is applied separately to 1-tag and > 2-tag 
samples and considers many different combinations of 



possible requirements. The smallest values for the 
uncertainty are obtained using (iV out > 0.90, x 2 < 6) 
in the 1-tag sample and (7V ut > 0.88, x 2 < 5) in the 
> 2-tag sample so that we add these requirements 
to the prerequisites described in Sec. IIII B[ The fi- 
nal definition of the samples used in our analysis is 
summarized in Table [TTJ 

After these selections, 3452 and 441 events are ob- 
served for the 1-tag and > 2-tag samples respec- 
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FIG. 4: Histograms of the reconstructed W mass m^ c for 
1-tag events, upper plot, and > 2-tag events, lower plot, 
are shown in a control region defined by 0.75 < 7V out < 
0.85. Along with the data are plotted the expected back- 
ground and the signal contribution for M top = 175 GeV/c 2 
and the default value AJES = 0, normalized to the differ- 
ence between the data and the background. The value of 
the purely statistical \ 2 probability is reported on each 
plot. 



tively. We can evaluate the average number of back- 
ground events expected in the selected samples and 
their uncertainties, as described in Sec. [V] The sys- 
tematic uncertainties on the background normaliza- 
tions are estimated by assuming that the discrepancy 
between the observed number of events in the data 
and the sum of the expected contributions from signal 
and background (where, in this case, the theoretical 
cross section value of 6.7 pb is considered for tt events 
production) is due to a bad evaluation of the back- 
ground. This is done separately for 1-tag and > 2-tag 
samples, and the resulting relative uncertainties on 
the expected number of events are o-(n, b ta ^,) = 2.9% 

and c^^cxpp = 14-6%, respectively. The efficien- 
cies of the full selection on tt events corresponding 
to M top = 175GeV/c 2 and AJES = are 3.6% and 
1.0% for 1-tag and > 2-tag events respectively. These 
values are used to evaluate the expected signal con- 
tributions of Table HIIl where o~ t i = 6.7 pb is assumed. 
In the same table, the observed number of events and 
the expected background in each sample are also sum- 
marized. 



TABLE III: Number of events observed in the selected 
data samples and corresponding expected numbers of 
background and tt events. The signal contribution is eval- 
uated for Mtop = 175 GeV/c 2 , AJES = and <j ti = 6.7pb. 



Event sample 


Observed 


Background 


tt 


One tag 


3452 


2785 ± 83 


693 


> 2 tags 


441 


201 ± 29 


193 



VIII. LIKELIHOOD FIT 

The technique described in Sec. lVIl allows one to 
obtain sets of observed m\ ec and m^>° values recon- 
structed in the data samples with 1 or > 2 tags as well 
as to build signal and background distributions for the 
same variables. In order to measure the top quark 
mass simultaneously with the JES, a fit is performed 
where an unbinned likelihood function is maximized to 
find the values of M top , AJES, and the number of sig- 
nal (n s ) and background (rib) events for each tagging 
category which give the probability density functions 
(p.d.f.'s) best describing the data. 



TABLE II: Final definition and requirements for selected 
event samples. 



Event Sample 


6-tags 


AUt 


mj ec fit x 2 


One tag 


= 1 


> 0.90 


< 6 


> 2-tags 


2 or 3 


> 0.88 


< 5 



A. Probability density functions 

The signal templates are fitted by normalized com- 
binations of Gamma and Gaussian p.d.f.'s, and the 
dependence of the shape on input Mt op and AJES is 
included writing the parameters of the p.d.f.'s as lin- 
ear functions of these variables. Figures [5] and [5] show 
examples of the fitted p.d.f.'s superimposed on the 
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FIG. 5: Histograms and corresponding fitted probability 
density functions for the signal m\ ec in > 2-tag events for 
a constant AJES value (AJES = 0) , varying the input top 
quark mass (upper plot) and for a constant Af top value 
(175GeV/c 2 ), varying the input jet energy scale (lower 
plot). 



FIG. 6: Histograms and corresponding fitted probability 
density functions for the signal m^ c in > 2-tag events for 
a constant AJES value (AJES = 0), varying the input 
top quark mass (upper plot), where the independence of 
m^ c on M top is apparent, and for a constant M top value 



(175 GeV/cr 
plot). 



), varying the input jet energy scale (lower 



m™ and rnyy signal templates respectively for differ- 
ent M top and AJES values. 

The shape of distributions built for the background 
cannot depend on the characteristics of signal events, 
and in particular on the value of top quark mass. 
Moreover, as they are obtained from data, the shapes 
correspond to the reference value of the jet energy 
scale. For these reasons no dependence on M top and 
JES is considered in the p.d.f.'s used to fit the back- 
ground templates. Actually, a very weak dependence 
is introduced through the corrections to the shape of 
the background distributions, performed to take into 



account the presence of signal events in the pretag 
sample, as described in Sec. lVI Cl These effects are 
taken into account as a systematic uncertainty. Ex- 
amples of background m\ cc and m^> c distributions and 
the corresponding fitted p.d.f.'s are shown in Fig. [7] 
for > 2-tag events. Discrepancies between the fitted 
p.d.f.'s and the corresponding distributions are consid- 
ered in the calibration procedure, presented in Sec. lIXI 
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FIG. 7: Data-driven background histograms of mj ec (up- 
per plot) and m^ c (lower plot) for > 2-tag events. The 
bands denote the la uncertainty on the bin contents of 
the histograms, including both statistical and systematic 
contributions. The solid lines show the p.d.f.'s fitted to 
the histograms. 



B. The Likelihood Function 



where the four terms on the right side assume respec- 
tively the following form [the superscripts referring to 
the tag sample are omitted and f s = n s /(n s + rib), 
fb = 1 - /,] : 

JVobB 

Hf*-P™g ( m t, i I Mtop, AJES) + f b - P££* (ith, i) , (14) 
i=l 

N obs 

I] /• • P Z W ( m w, * I M to P , AJES) + h ■ P™ k g (™w, i) , (15) 

(16) 



e -(n.+n b ) . , + m) N ob 



exp 



iVobs! 



{n b — n( ftj exp))' 



(17) 



In expression (|14|) the probability to observe the set 
mt,i, (i = 1,..., iVobs) of mj cc values reconstructed 
in the data is calculated by using the total probabil- 
ity density function resulting from the combination 
of the parametrized signal and background p.d.f.'s 
(Sec lVIIlX]) . and P^* respectively, as a func- 

tion of the free parameters of the fit. In term ([T5]) the 
same is done for the set of the observed W masses, 
fnw,h (i = 1, •••) Aobs), and the m^ c p.d.f. The 
term (IT5|) . £p i SS , gives the probability to observe the 
number of events selected in the data, given the aver- 
age number of signal (n s ) and background (rib) events 
expected in the sample, as assumed at each step of 
the likelihood fit. In the last term, (fT7|) . the param- 
eter rib is constrained by a Gaussian to the a priori 
background estimate given in Sec. IVII1 i.e. cxp ) = 
2785 ± 83 for 1-tag events and ntb,exp) = 201 ± 29 
for > 2-tag events. Finally, the last term in expres- 
sion (jT2"j) , Ca jES C o„ s tr * is a Gaussian term constraining 
AJES to its a priori value : 



exp 



(AJES - AJES 

const r ) 



(18) 



When the measurement is performed on data, the JES 
can be constrained to the value independently mea- 
sured in 16]. Given the meaning of AJES, described 
in Sec. lIVl this means that, in this case, AJES cons t r = 
0. 



The likelihood function C is divided into three main 
parts and can be written as 

£ = £l tag X £>2 tags X ^AJES oonstr • (12) 

The L\ and £>2tags terms further consist of other 
factors : 

£l,>2tags = -CjUtop x -^JES x ^Poiss X /^bkg , (13) 



IX. VERIFICATION AND CALIBRATION OF 
THE METHOD 

We want to investigate the possible presence of bi- 
ases in the top quark mass and jet energy scale mea- 
surements introduced by our method, as well as to 
have an estimate of its statistical power before per- 
forming the measurement on the actual data sam- 
ple. To do so, we run realistic PEs where pseudo- 
data are extracted from simulated signal and data- 
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driven background distributions. A set of 3000 PEs is 
performed for each simulated value of the top quark 
mass and of the displacement in the jet energy scale 
(Sec. lIVp . Using the notation introduced in SeclIV) we 
refer to these input values as M t '" p and AJES m and 
they represent the true values we want to measure. In 
each PE the actual numbers of signal (-/Vr S)0 & s )) and 
background (N^ b s ) ) events in each tagging category 
are generated with Poisson distributions with mean 
«( s ,cx P ) = - oxp ) and n (bi0xp) respectively, 

where iV bs a re the observed number of events in the 
data samples (A b s = 3452 for 1-tag and A b s = 441 
for > 2-tag events). A set of A( Si0&s ) and A(b !0&s ) 
mass values is then drawn from m\ cc and m 1 ^ distri- 
butions of signal and ba ckgrou nd and used as input 
to the likelihood fit (Sec. lVIlfl) that returns simulta- 
neous measurements of M top and AJES, denoted as 
Af t ° U p and AJES out . The average of these measure- 
ments over the whole set of 3000 PEs represents the 
best estimate of the input values obtained by the fit- 
ting procedure and therefore can be used to study its 
behavior. We fit the dependence of these averages 
with respect to the input values over the whole range 
of simulated M™ p and AJES in as 

<AO = (A 00 + Aoi ■ AJES in ) 

+ (A 10 + An ■ AJES in ) • (AC P - 175) , (19) 



TABLE IV: Coefficients of expressions HU) and ((20j) as 
obtained from calibrated and uncalibrated pseudoexperi- 



ments. 


The ideal values in 


the absence of any bias 


are also 


shown. 










Uncalibrated 


Oali bratecl 


iucai 




PEs 


PEs 


value 


Ann 
J -HJU 


175 47 ± 01 


174.99 ±0.01 


175 


A i 


-0.24 ±0.01 


0.00 ±0.01 





Aio 


0.985 ±0.002 


1.000 ±0.002 


1 


An 


0.009 ±0.001 


0.000 ±0.001 





-Boo 


-0.026 ±0.003 


0.002 ±0.003 





-Boi 


0.0009 ±0.0004 


-0.0001 ±0.0004 





-Bio 


1.052 ±0.002 


0.999 ±0.002 


1 


B n 


-0.0016 ±0.0002 


0.0001 ± 0.0002 






sured top quark mass and jet energy scale displace- 
ment for true M top and AJES around 175 GeV/c 2 and 
are 

(5M t c ° p rr (stat + JES) = 2.0GcV/c 2 , (21) 
(5AJES corr (stat + M top ) = 0.45. (22) 



(AJES out ) = [Boo + B 01 ■ (AC P - 175)] 

+ [b u + Bn ■ (AC P ~ 175)] • AJES in . (20) 

These relations can be inverted to obtain calibration 
functions to be applied to further measurements and 
therefore, on average, a more reliable estimate of the 
true values (2D calibration) . The calibrated values re- 
sulting from a measurement giving M£Hj and AJES out 
are denoted as Af t c ° p r and AJES corr , while the respec- 
tive uncertainties, obtained by propagating through 
the calibration the uncertainties from the likelihood 
fit, are (5M t c ° p rr and <5AJES corr . A second set of PEs is 
then performed to test the goodness of the procedure. 
Table IIVI shows the coefficients Ay and obtained 
both from calibrated and uncalibrated PEs compared 
to their ideal values in the absence of any bias. 

In Fig. [5] examples of linearity plots are shown for 
calibrated PEs. These plots, together with the num- 
bers in Table IIV( show how the calibration removes 
any average bias. To check that the uncertainties 
5M£°™ and (5AJES corr are also unbiased, we consider 
the width of Mjr™ 1 and AJES corr pull distributions, 
that is, the distributions of deviations of the calibrated 
values from the true inputs in the PEs, divided by 
the uncertainties themselves. We find that the uncer- 
tainties are both underestimated, and multiplicative 
correction factors equal to 1.084 for (5M t c ° p r and to 
1.115 for <5AJES corr are needed. After these correc- 
tions the average expected uncertainty on the mea- 



X. SYSTEMATIC UNCERTAINTIES 

The possible systematic uncertainties on the top 
quark mass and the jet energy scale measurements 
have been studied and are summarized in this sec- 
tion. These arise mostly from the measurement tech- 
nique itself, from uncertainties in the simulation of 
the ti events, from mismodeling of the detector re- 
sponse, and from uncertainty on the shapes of signal 
and background templates used to derive the p.d.f.'s 
and to calibrate the measurement. The uncertainties 
are usually evaluated by performing PEs, extracting 
pseudodata from templates built using signal samples 
where the possible systematic effects have been con- 
sidered and included. Corrections to the shape of the 
raw background templates are performed as described 
in Sec. lVICl to obtain the corrected background tem- 
plates corresponding to the effect one wants to study. 
On the contrary, nothing is changed in the elements of 
the likelihood fit, because it is the default procedure 
that we want to apply to real data and that, there- 
fore, we have to test in case of possible mismodeling 
of the data themselves. The results from these PEs 
are then compared to the ones obtained by using de- 
fault templates, and the shifts in the average Af t c ° p r 
and AJES corr values are taken as the estimate of the 
systematic uncertainties. In some cases the statisti- 
cal uncertainties on the shifts may be larger than the 
shifts themselves, and therefore we use conservatively 
the former as the systematic uncertainty. In the fol- 
lowing, after the description of each effect, we also 
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FIG. 8: Examples of response linearity plots ((Af t C op r ) 
vs M£ p , upper plot, and <AJES corr ) vs AJES in , lower 
plot) after the 2D calibration. The solid lines represent the 
linear functions which best fit the response as a function of 
the input values, while the dashed straight lines represent 
the ideal behavior. 



quote in parentheses the values of the corresponding 
uncertainties for the top quark mass and the jet energy 
scale respectively. These values are then summarized 
in TableEl 

The 2D calibration removes the average biases, es- 
pecially related to the parametrization of the tem- 
plates using smooth probability density functions. 
Residual biases usually exist at single (M"^ p , AJES m ) 
points, and have to be taken into account. We there- 
fore consider the shift of the mean of the pull distribu- 
tions with respect to at each (M t '£ , AJES 111 ) point 
to evaluate this residual bias systematic uncertainty, 



which, given the definition of pull in Sec. lIXl is defined 
as a function of the uncertainty on the calibrated mea- 
surements. To obtain the proper coverage of both pos- 
itive and negative biases we evaluate them separately, 
so that asymmetric uncertainties are finally consid- 
ered. They are generally given by (lollo) ' ^-^top 1 f° r 
M™" and (t°;^)-<5AJES corr for AJES corr . Specifying 
the values obtained in the measurement on the data, 
described in Sec. ED we obtain GeV/c 2 on M t c ° p rr , 

t° H on AJES corr . 

The uncertainties on the parameters of the 2D cali- 
bration give a small uncertainty on the corrected val- 
ues M™* T and AJES corr which can be evaluated by 
the calibration functions and the values of Af top and 
AJES fitted in the data (< 0.1 GeV/c 2 , < 0.01). 

Many sources of systematic effects arise from un- 
certainties in modeling of the hard interaction in sim- 
ulated events, pythia and HERWIG [11] Monte Carlo 
generators differ in their hadronization schemes and in 
their description of the underlying event and multiple 
interactions. The default signal samples have been 
generated with PYTHIA and therefore an uncertainty 
is obtained by using a sample generated using herwig 
(0.3GeV/c 2 , 0.25). 

Jets coming from possible emission of hard gluons 
might fall among the six leading jets and populate the 
tails in the top quark invariant mass distribution. The 
amount of radiation from partons in the initial or final 
state is set by parameters of the pythia generator 
used to simulate signal events. To study these effects, 
templates are built using samples where the values 
of the parameters have been changed with respect to 
the default, to increase or to decrease the amount of 
radiation Q (0.1GeV/c 2 , 0.06). 

Since the default jet energy corrections are derived 
from data samples deficient in heavy flavors an 
additional uncertainty comes from considering the dif- 
ferent properties of b quarks. We account for the 
uncertainties on the 6-quark semileptonic branching 
ratios, fragmentation modeling, and calorimeter re- 
sponse to heavy-flavor hadrons (0.2GeV/c 2 , 0.04). 

The different efficiency of the 6-tagging algorithm 
on data and simulated events is usually considered by 
introducing a constant scale factor (6-tag SF). The 
overall uncertainty on this parameter affects the cross 
section measurement described in Sec. lXIII However, 
such scale factor does not need to be considered re- 
garding the top quark mass measurement, because it 
could slightly change only the population of the signal 
templates, but not their shape. On the other hand, 
variations of the latter could be caused by the possible 
dependence of the 6-tag SF on the transverse energy 
of jets, which is then considered as a systematic effect 
(0.1GeV/c 2 , 0.01). 

The uncertainty on the top quark mass coming from 
the likelihood fit includes the uncertainty due to the 
jet energy scale. However, as described in Sec. lIII Bl 
this uncertainty is the result of many independent ef- 
fects with different behavior with respect to proper- 
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ties of jets like Et and r\ |16|, and therefore represents 
a leading-order correction. Second-order effects can 
arise from uncertainties on the single corrections ap- 
plied to the jet energies. To evaluate these possible ef- 
fects, we build signal templates by varying separately 
by ±1(7 the single corrections and, for each one of 
these variations, PEs were performed by using these 
templates and not applying the constraint £AjES C o„ 3 tr 
in the likelihood fit, as this term is related to effects 
of the full correction. The resulting uncertainties have 
been added in quadrature to obtain a residual JES un- 
certainty on the top quark mass (0.5GeV/c 2 ). 

The choice of parton distribution functions (PDF) 
in the proton used in Monte Carlo generators can af- 
fect the kinematics of simulated tt events and thus 
the top quark mass measurement. We considered 
four sources of uncertainties : the difference arising 
from the use of the default CTEQ5L H3 PDF and 
one calculated from the MRST group, MRST72 [H]; 
the uncertainty depending on the value of a s , eval- 
uated by the difference between the use of MRST72 
and MRST75 PDF's; the uncertainty depending on 
the differences between the LO and NLO calculations 
of PDF' s, evaluated by the difference between using 
default CTEQ5L (LO) and CTEQ6M (NLO) PDF's; 
and the uncertainties on PDF' s derived from experi- 
mental data uncertainties (l~ ; 2 GeV/c 2 , to.ol)- 

The probability to have multiple pp interactions 
during the same bunch crossing is a function of the 
instantaneous luminosity. This is reflected in the in- 
creasing number of primary vertices reconstructed in 
the events at higher luminosities. We account for the 
fact that the simulated samples for the signal pro- 
cess do not model the actual luminosity profile of the 
data by considering the signal distributions for events 
with 1,2,3, and > 4 reconstructed vertices separately. 
These distributions are then used to obtain the tem- 
plates by weighted averages, where the weights are 
evaluated as the fractions of events with 1, 2, 3 and 
> 4 vertices observed in the data. Moreover, a possi- 
ble mismodeling of the dependence of the jet energy 
response as a function of the reconstructed number 
of primary vertices in simulated events is considered 
(0.2GeV/c 2 , 0.01). 

Uncertainties from modeling of color reconnection 
effects are estimated by comparing the results of 
two sets of PEs performed drawing pseudodata from 
templates built using two different samples of events 
simulated by PYTHIA. The samples are generated with 
two different sets of parameters, corresponding to two 
different models of color reconnection (0.4GeV/c 2 , 
0.08). 

The shapes of the signal and background distribu- 
tions are affected by uncertainties due to the limited 
statistics of the simulated events and data samples 
used to build them. These uncertainties affect the 
results of a measurement which is performed max- 
imizing an unbinned likelihood where parametrized 
p.d.f.' s, fitted to default templates, are evaluated. We 



address this effect obtaining 100 sets of templates by 
statistical fluctuations of default ones and performing 
pseudoexperiments drawing data from each of these 
sets separately. From each set we obtain an aver- 
age value for M t c " r and AJES corr , and the spread of 
these values is taken as the systematic uncertainty 
(0.3 GeV/c 2 , 0.07). 

Besides the purely statistical effects, quoted above, 
the shape of the background templates also has un- 
certainties due to the corrections for the presence of 
signal events in the pretag sample, Sec. IVICl and to 
the systematic uncertainty on the background normal- 
ization, Sec. I VIII We address this source of system- 
atic uncertainty by the same technique used for the 
statistical contributions, that is by obtaining a set of 
100 background templates where the content of each 
bin is separately fluctuated by Gaussian distributions 
centered on the default bin content and with a width 
equal to its systematic uncertainty, and taking the 
spread of results from PEs as the systematic uncer- 
tainty (0.1GeV/c 2 , 0.02). 

Table [V] shows a summary of all the systematic un- 
certainties and their sum in quadrature, which gives 
a total systematic uncertainty of t-ix,GeV/c 2 for the 
Mtop measurement and Iq^* for the AJES. 

TABLE V: Systematic uncertainties and their sizes for 
the top-quark mass and the jet energy scale measurements. 
The total uncertainty is obtained by the sum in quadrature 
of single contributions. 



Source 


5M£? 
(GeV/c 2 ) 


SAJES syst 


Residual bias 


+0.8 

-0.4 


+0.18 
-0.24 


2D calibration 


< 0.1 


< 0.01 


Generator 


0.3 


0.25 


ISR/FSR 


0.1 


0.06 


6-jet energy scale 


0.2 


0.04 


6-tag SF E T dependence 


0.1 


0.01 


Residual JES 


0.5 




PDF 


+0.3 
-0.2 


+0.05 
-0.04 


Multiple pp interactions 


0.2 


0.01 


Color reconnection 


0.4 


0.08 


Statistics of templates 


0.3 


0.07 


Background shape 


0.1 


0.02 


Total 


+1.2 
-1.0 


+0.34 
-0.37 



XI. TOP MASS AND JES MEASUREMENTS 

After the kinematical selection with N ou t > 0.90 
(> 0.88) and x 2 < 6 (< 5), we are left with 3452 
(441) events with one (> 2) tag(s). The background 
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amounts to 2785 ± 83 (201 ± 29) for events with one 
(> 2) tag(s). 

For these events a top quark mass has been 
reconstructed using the likelihood fit described in 
Sec. IVIII Bl and applied to the data sample. Once 
the calibration procedure and corrections are applied, 
as described in Sec. IIX1 the best estimate of the top 
quark mass is 

M top = 174.8 ± 2.4 (stat + JES) GeV/c 2 , (23) 

while the value obtained for the jet energy scale dis- 
placement is 



AJES = -0.30 ± 0.47 (stat + M top ) . 



(24) 



We can also evaluate separately the purely statistical 
contributions obtaining 

Mtop = 174.8 ± 1.7 (stat) ± 1.6 (JES) GeV/c 2 , (25) 

and 

AJES = -0.30 ± 0.35 (stat) ± 0.32 (M top ) . (26) 

The plot in Fig. [9] shows the measured values to- 
gether with the log-likelihood contours corresponding 
to 1 a, 2 a, and 3 a uncertainty on the value of the top 
quark mass[24j]. The slope of the major axis of the 
contours denotes that the measurements of M top and 
AJES have a negative correlation, and the value of the 
correlation coefficient obtained from the likelihood fit 
0.68. 



is 



The plots in Fig. [TU] show the mji cc distributions 
for the data compared to the expected background 
and the signal for a top quark mass of 175.0 GeV/c 2 
and a jet energy scale displacement of —0.3, that is, 
the values of simulated M top and AJES as close as 
possible to the measurements in the data. The signal 
and background contributions are normalized to the 
respective number of events as fitted in the data. 

The plots in Fig. [11] compare the measured statis- 
tical uncertainty, just after the 2D calibration, with 
the expected distribution from default pseudoexper- 
iments using as inputs A/top = 175.0 GeV/c 2 and 
AJES = —0.3. We find that the probability of achiev- 
ing a better sensitivity is 91.6% for M top and 81.2% 
for AJES. 

Summarizing, the measured values for the top quark 
mass and the jet energy scale are 
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FIG. 9: Negative log-likelihood contours for the likeli- 
hood fit performed for the A/ top and AJES measurements. 
The minimum is also shown and corresponds to the values 
measured in the data. The contours are drawn at values of 
0.5, 2.0, and 4.5 of the increase of the log-likelihood from 
the minimum value. These curves correspond to la, 2a, 
and 3 a uncertainty on the measurement of the top quark 



written as 



M top = 174.8 ± 1.7 (stat) l(% (syst) GeV/c 2 , (29) 



AJES = -0.30 ± 0.35 (stat) toil ( s Y st ) • ( 30 ) 

This measurement of the top quark mass has been 
used in the current world average |5[ . 



XII. CROSS SECTION MEASUREMENT 

The procedure used to measure the top quark mass 
also returns the average number of signal events ex- 
pected, given the selected data samples. These results 
can be turned into a measurement of the tt cross sec- 
tion, as follows. 



The likelihood function 



M 



top 



174.8 ± 2.4 (stat + JES) Z l (syst) GeV/c 2 , (27) 



AJES = -0.30 ± 0.47 (stat + M top ) (syst) , (28) 

which, isolating the purely statistical contributions 
and adding the uncertainties from JES and M top to 
the respective systematic uncertainties, can also be 



From the number of signal events, n l s tag and 
n^ 2tags , as obtained from the mass likelihood fit, we 
derive a measurement of the tt production cross sec- 
tion considering the efficiency for selecting a tt event 
in the two tagging categories. 

The cross section measurement is performed by 
maximizing a likelihood function which can be divided 
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FIG. 10: Histograms of ml ec as obtained in the data (black 
points) for 1-tag (upper plot) and > 2-tag events (lower 
plot) are compared to the distributions from signal and 
background corresponding to M top = 175 GeV/c 2 and 
AJES = —0.3. The expected histograms are normalized 
to the measured values for the average number of signal 
and background events. The values of the purely statis- 
tical x 2 an d of its probability are reported on each plot, 
showing the overall agreement between the data and the 
distributions corresponding to the fitted values of M top 
and AJES. 



into two parts : 

£ = £ 

1 tag * £>2tagsi (31) 

where each term can be expressed as : 

£l,>2tags = £<j tl X £e > (32) 



M fnn = 175.0 GeV/c 2 AJES = -0.3 

top 



Expected 
Measured 




1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 3.2 

5M, (stat + JES) [GeV/c 2 ] 



CM F 
<=> 0.3- 
© 
W 
HI 

0- 0.25 - 



U 
(0 



M, = 175.0 GeV/c 2 AJES = -0.3 



Expected 



Measured 



0.15- 



0.1- 



0.05- 




. I i i . . I . . i . I . i . . I 



I . i . . I . . i i I . . . . I . 



0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 

5AJES(stat + M ) 

v top' 



FIG. 11: Distributions of the uncertainties on the top 
quark mass (upper plot) and the jet energy scale dis- 
placement (lower plot) as expected from default PEs per- 
formed using as input values Mto P = 175.0 GeV/c 2 and 
AJES m = —0.3. The vertical lines indicate the uncertain- 
ties obtained in our reported result. 



where 



exp 



(a t i- e ■ L - n s y 
2ai 



(33) 



contains all the parameters of the fit, i.e. the produc- 
tion cross section cr t j , the integrated luminosity L, the 
signal efficiency e, and the signal yield n s ± a ns , as 
given by the mass measurement, while C e is a Gaus- 
sian term constraining the efficiency within its statis- 
tical uncertainty. 

The efficiencies are evaluated using a sam- 
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TABLE VI: Input variables to the cross section evaluation. 
For the signal yields, the first uncertainty is the purely 
statistical one. 



Variable 



Input value 



Signal yield, one tag 
Signal yield, > 2 tags 
Efficiency, one tag 
Efficiency, > 2 tags 
Integrated luminosity 



643 ± 59 ± 54 
216 ±21 ±14 
(3.55 ±0.01)% 
(1.00 ±0.01)% 
2874 ± 172 pb- 1 



TABLE VII: Cross section as evaluated assuming different 
values for M top and AJES. 



M top (GeV/c 2 ) 


AJES 


&tt (P b ) 


175.0 


-0.3 


7.24 


175.0 


0.0 


7.00 


172.5 


0.0 


7.21 


170.0 


0.0 


7.29 



pie of about 4 x 10 6 tt events generated with 
M top = 175GeV/c 2 and assuming AJES = -0.3, i.e. 
the value we measured by the mass likelihood fit, and 
are summarized along with signal yields and other pa- 
rameters in TablelVTl 



While studying the performance of the procedure, 
using pseudoexperiments produced assuming a given 
input cross section, we observe the need to introduce 
a small correction. The outcome of the fit needs to be 
multiplied by a factor k a — 0.982 ± 0.008 in order to 
obtain an unbiased measurement of the cross section. 

From the maximization of the likelihood, we obtain 
a central value for the tt production cross section 



o tl = 7.2 ± 0.5(stat) ± 0.4(lum) pb , 



(34) 



evaluated assuming Mtop = 175 GeV/c 2 and AJES = 
—0.3, close to the values measured in Sec. lXIl The 
first uncertainty is the statistical one, while the sec- 
ond one derives from the 6% uncertainty on the in- 
tegrated luminosity. As the signal efficiencies depend 
strongly on the assumed values for M top and AJES, 
the measured tt cross section also has the same depen- 
dence. For reference we report in Table lVlIl the cross 
sections corresponding to other (M top ,AJES) points 
with a top quark mass near the current CDF average. 
In this case we assume AJES = and the systematic 
uncertainty on JES is increased from 6.1% to 9.2%, 
corresponding to changing the AJES by ±1 rather 
than by ±0.6 units, that is the sum in quadrature 
of the uncertainties on the measured jet energy scale, 
SecKU 



TABLE VIII: Systematic uncertainties and their relative 
sizes for the cross section measurement. The total un- 
certainty is obtained by the sum in quadrature of single 
contributions. 



Source 


Afftf/ow (%) 


Calibration factor 


0.8 


Generator 


4.2 


ISR/FSR 


0.6 


6-jet energy scale 


2.8 


6-tag SF Et dependence 


5.4 


PDF 


3.4 


Multiple pp interactions 


2.5 


Color reconnection 


0.8 


Templates statistics 


0.8 


Background shape 


0.3 


Background normalization 


8.2 


JES 


6.1 


Residual JES 


2.1 


Primary vertex 


0.2 


BR(W^ -> hadrons) 


0.8 


Trigger 


1.8 


Total 


13.7 



B. Systematic uncertainties 



Most of the sources of systematic uncertainties af- 
fecting the measurement of a t i are the same discussed 
for the measurement of the top quark mass. We just 
need to evaluate their effects both on the signal yields 
and on the signal efficiencies in order to derive the ef- 
fects on the cross section. There are few other sources 
of systematic uncertainty specific to a cross section 
measurement which have not been discussed in Sec. IXl 
because they affect only the signal efficiencies. These 
include the uncertainty on the calibration constant, 
k a , on the W — > hadrons branching ratio (BR) (24(, 
on the trigger simulation and on the distribution of 
the primary vertex z-coordinate. As for the effect of 
the JES uncertainty on the efficiency, we have eval- 
uated it by changing the AJES by ±0.6 units with 
respect to the measured value AJES = —0.3. Resid- 
ual effects due to individual levels of corrections have 
been accounted for, too. The relative uncertainties 
Aa t i/&tt for the individual sources are summarized in 
Table lVIlTl Considering their sum in quadrature, the 
tt production cross section amounts to 

cr tt -= 7.2 ± 0.5(stat) ± l.O(syst) ± 0.4(lum) pb , (35) 

assuming M top = 175 GeV/c 2 and AJES = -0.3. 
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XIII. CONCLUSIONS 

Using a very effective neural-network-based kinc- 
matical selection and a b jet identification technique, 
we measure the top quark mass to be 

M t0 p = 174.8 ± 2.4(stat + JES) ti'.o(syst) GeV/c 2 , (36) 

and the tt production cross section to be 

o t i = 7.2 ± 0.5(stat) ± l.O(syst) ± 0.4(lum) pb . (37) 

These values represent the most precise measure- 
ments to date of -Mt op and a t i in the all-hadronic decay 
channel. The results are consistent with the measure- 
ments obtained in other decay channels by CDF and 
DO Collaborations || 0| and, as it concerns cr t j, with 
the theoretical predictions evaluated at the value of 
the top quark mass obtained in our measurement Q . 
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